Stationary dark energy with a baryon-dominated era: solving the coincidence problem 

with a linear coupling 
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We show that all cosmological models with an accelerated stationary global attractor reduce 
asymptotically to a dark energy field with an exponential potential coupled linearly to a perfect 
fluid dark matter. In such models the abundance of the dark components reaches a stationary value 
and therefore the problem of their present coincidence is solved. The requirement of a vanishing 
coupling of the baryons in order to pass local gravity experiments induces the existence of an inter- 
mediate baryon-dominated era. We discuss in detail the properties of these models and show that 
to accomodate standard nucleosynthesis they cannot produce a microwave background consistent 
with observations. We conclude that, among stationary models, only a time-dependent coupling or 
equation of state might provide a realistic cosmology. 
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INTRODUCTION 



The recent observations of an accelerating expansion 
0, together with new cosmic microwave, large-scale 
structure and lensing data, give a strong indication that 
the universe fluid is composed of at least four different 
components: relativistic matter (<0.01%), baryons (few 
per cent), dark matter (~30%) and dark energy (~70%). 
The present amount of these components raises some 
deep questions: Why are the dark matter and dark en- 
ergy, which supposedly have a different scaling with time, 
almost equal right now Why are the baryons strongly 
suppressed with respect to dark matter? Why is the coin- 
cidence between the dark components relatively close to 
the equivalence between matter and radiation (as pointed 
out by §)? 

The problem of the present coincidence between dark 
energy and dark matter could be simply solved if the two 
fluids have the same scaling with time (let us call a sys- 
tem of fluids p\ , p2 with identical scaling a "stationary" 
model, since d{p\j p<i)j dt = 0). It is interesting to remark 
that such a scaling can soon be observationally tested |Q. 
It is well known that, assuming an exponential potential 
for the scalar field representing the dark energy, the field 
scales as the dominant component || |[ [t|] , but this does 
not produce acceleration. The question therefore arises 
of which is the most general system of dark matter (mod- 
eled as a perfect fluid) and dark energy (a scalar field) 
that allows an accelerated stationary global attractor. In 
this paper, following the arguments of ref. J^], we show 
that all models which contain an accelerated stationary 
global attractor reduce asymptotically to a system char- 
acterized by an exponential potential and a linear cou- 
pling between the two components. Systems of such a 
kind have been already analyzed first in refs. ||_1Q] and 
successively by many authors [|n], [l2], [l3[ [l^, J15 [while 
their perturbations have been studied in ref. []16] . In or- 
der to pass local gravity experiments jl7), we argue that 
it is necessary to break the universality of the coupling, 
leaving the baryons uncoupled or weakly coupled (see 



e.g. [I?], |T|, |(J). These assumptions completely 
define our cosmology. 

The consequences of a linear coupling between the dark 
components are manifold. First, as already remarked, 
this explains the cosmic coincidence and the acceleration. 
Second, the accelerated regime explains the decay of the 
baryons with respect to the dark components. Third, the 
near coincidence of the equivalence between the luminous 
components (baryons and radiation) and the beginning 
of the dark era (dark equivalence from now) is automat- 
ically enforced. Fourth, the baryons are the dominant 
component between the two equivalence times, produc- 
ing a decelerated epoch in which gravitational instability 
is effective. Fifth, the present status of the universe is 
independent of the initial conditions, being on a global 
attractor. In addition, it is to be noticed that the model 
requires only constants of order unity in Planck units, 
and that they are all fixed by the observations of the 
present energy densities and the acceleration. 

Nothwithstanding these intriguing features, we show 
that the model we present here fails in satisfying at 
the same time the nucleosynthesis requirements together 
with producing an acceptable cosmic microwave back- 
ground (CMB) angular power spectrum. In fact if the 
conditions for a standard nucleosynthesis are adopted, 
during the recent accelerated regime a fast growth of the 
perturbations is induced, which in turn causes an ex- 
cessive integrated Sachs- Wolfe (ISW) effect, in contrast 
with the observations. Nevertheless, we think that the 
dynamics we discuss is interesting on its own. 

The conclusion we draw is that only models with a 
time-dependent coupling or a potential more complicated 
that an exponential may contain an accelerated global 
attractor and at the same time be compatible with nu- 
cleosynthesis. A solution along these lines has been pro- 
posed in j2l| where a non-linear modulation of the cou- 
pling allows nucleosynthesis to happen and structure to 
form in a regime of weak coupling, while the acceleration 
is produced in the subsequent (present) regime of strong 
coupling. In ref. it has been shown that a similar 
mechanism may be realized in superstring theories. 
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II. HOMOGENEOUS SOLUTIONS 



with 



Let us show first why all two-fluid systems with a sta- 
tionary accelerated attractor reduce asymptotically to 
the one investigated below. Let us write a generic coupled 
two-fluid systems with equations of state p x — (w x — l)p x 
in a flat Friedmann metric as 



Pc 



- 3Hw c p c 
3Hw^p^ 



8. 



-s, 



(1) 

(2) 



where the subscript c stands for cold dark matter (CDM) 
and the subscript 4> f° r a scalar field. The Friedmann 
equation is 

3H 2 = k 2 (pj, + p c ) , 

where k 2 — 8ir and G = c = 1. As shown in ref. the 
stationary condition d(p c / p$)dt — 0, that is = Ap c , 
can be satisfied only if 



A 



w<t>) 



fPc 



(3) 



Putting w c 
have 



1 and observing that A = 0^/(1 — O^), we 



3ft 2 i!^ 77 — 1 
2 VTTv 



\<j>\Pc 



(4) 



where 77 = 2U/(j) 2 , the ratio of the potential to kinetic 
scalar field energy (notice that — 2/(1 + 77)). The 
stationary solution is accelerated if 



a - = H- 
a 



H 2 = k 2 (1 + A)p c 



2 6 



>0, 



which, for 2 > > 0, can be realized only if fl$ > 1/3 
(i.e. A > 1/2) and 



r? > 



3ft,* - 1 



> 2. 



(5) 



So far we repeated the steps of ref. ||. Now, let us 
consider the asymptotic behavior of 77. If 77 — * 0, the 
kinetic energy dominates over the potential energy, and 
the asymptotic solution is not accelerated. If, on the 
other hand, 77 — > 00, the potential energy dominates, and 
the solution will be identical to that of a cosmological 
constant. In fact, in this limit qb — and w^, — and Eq. 
(0) gives 8 — where 



8 = p c J3 K 2 U((/))n^ 



(6) 



which implies p c — » 0. Therefore, barring oscillatory so- 
lutions, only if 77 — > const, the scalar field behaves as 
stationary accelerated dark energy (clearly this case in- 
cludes also that of dark energy as a perfect fluid) . When 
77 is constant, the coupling reduces to 



8 = ^2/3k{3\<P\ Pc 



(7) 



(«) 



which is the form we study below. Moreover, it is easy 
to show that 77 =const. implies an exponential potential 
U = Uoe-V^^ where 



/' 



l + i(r7-l)(l-O ) 



(9) 



The conclusion is that a linear coupling and an exponen- 
tial potential represents the only non-trivial asymptotic 
case of stationary accelerated dark energy. It is not dif- 
ficult to see that this theorem extends also to a Brans- 
Dicke theory with an explicit coupling between matter 
components. As will be shown below, such a solution is 
also a global attractor in a certain region of the parame- 
ter space. 

The condition for acceleration, 77 > (30,^+1) / (30,^— 1), 
together with f2^ < 1, imply in Eq. (||) that 



(3 > V3/2. 



(10) 



This limit is much larger than allowed by local gravity 
experiments on baryons, which give at most < 0.01 
(see e.g. |^, [l7)). Therefore, the theory must break the 
universality of the coupling and let the baryons be de- 
coupled from dark energy. An immediate consequence 
of the species-dependent coupling, so far unnoticed, can 
be seen by observing that the energy density of the dark 
components scales as 



P4> ~ Pc 



(11) 



For any (3 > the energy density decays slower than 
in the standard matter-dominated Friedmann universe. 
Therefore, any uncoupled (or weakly coupled) compo- 
nent, as the baryons, decays faster than the coupled ones 
(see also @). 

The cosmology we study below is a more realistic ver- 
sion of the one above: we include in fact radiation and 
baryons, both of which are coupled to the dark compo- 
nents only through gravitation. Once baryons and radi- 
ation decay away, we recover the stationary accelerated 
attractor. The Einstein equations for our model have 
been already described in p| , in which a similar model 
(but on a different attractor, i.e., for different parame- 
ters) was studied (see also ref. |2^|). Here we summarize 
their properties. The conservation equations for the field 
4> , cold dark matter, baryons (b ), and radiation (7 ), 
plus the Friedmann equation, are 



- 3H</> + 

p c + 3Hp c 
p b + 3Hp b 
p y + 4iJp 7 



3H 2 = 




(12) 
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where H = a/a and U(<p) = [/ oe -V 2 / 3 ^ ( we put ^ in . 
stead of |</>| for generality). The coupling (3 can be seen 
as the relative strength of the dark matter-dark energy 
interaction with respect to the gravitational force. The 
only parameters of our model are (3 and li (the constant 
Uq can always be rescaled away by a redefinition of </>). 
For (3 = n = we reduce to the standard cosmological 
constant case, while for /3 = we recover the Ferreira & 
Joyce model of ||. As shown in ref. [Q, the coupling we 
assume here can be derived by a conformal transforma- 
tion of a Brans-Dicke model, which automatically leaves 
the radiation uncoupled. To decouple the baryons one 
needs to consider a two-metric Brans-Dicke Lagrangian 
as proposed in [jl7j. Additional theoretical motivations 
for this kind of coupling have been put forward in ref. 
ifHl and for coupled dark energy in general in ref. | ^ |. 
The system ( |12| ) is best studied in the new variables 

§y/U/3, z=§ ^p~/3 and 



y 



u = and the time variable a — log a. Then we 

obtain 

x = (z j z — 1) x — jiy 2 + (3(1 — x 2 — y 2 — z 2 — u 2 ), 

y' = pxy + y(2 + z'/z), 

v! = -3/2u + u(2 + z'/z) , (13) 



-z 1-31' 



zy 2 



A 



where the prime denotes derivation with respect to a. 
The CDM energy density parameter is obviously Vl c — 
1 — x 2 — y 2 — z 2 — u 2 while we also have 0,$ = x 2 + y 2 , 



Cl~ = z 2 and fib = 



The system is subject to the 



condition x 2 + y 2 + z 2 + u 2 < 1. 

The critical points of system ( fl^ ) are listed in Tab. I, 
where p is the scale factor exponent, a ~ T vl^-p = jp j 
where g = 4/3 2 + 4/3 n + 18, and where we used the sub- 
scripts b, c, r to denote the existence of baryons, matter 
or radiation, respectively, beside dark energy. In Tab. II 
we report the conditions of existence and stability of the 
critical points, denoting /i+ = (—0 + y/18 + /3 2 )/2 and 

Mo = -P- 20- 

In Fig. 1 we display the parameter space of the model, 
indicating for any choice of the parameters which point 
is a global attractor (notice that there is complete sym- 
metry under p — ► — /3 and [i — > — /x). As in ||ll[, in which 
the baryons have been included only as a perturbation, 
there exists one and only one global attractor for any 
choice of the parameters. The explicit inclusion of the 
baryons induces here two new critical points (bb and fb 
); moreover, contrary to [0, all the critical points with 
non-vanishing radiation are always unstable. 

From now on, we focus our attention on those param- 
eters for which the global attractor is b c , the only critical 
point that may be stationary and accelerated. In Fig. 1 
we show as a grey region the parameters for which this 
attractor is accelerated. When b c is the global attractor, 
the system goes through three phases: 

a) the radiation dominated era (the saddle b r ); 

6) the baryon dominated era ( the saddle bb ); 



c) the dark energy era (the global attractor b c ). 

The dynamics of the model is represented in Fig. 2 
(trend of £l c ,b,i,<t>) and in Fig. 3 (w e ff). During the 
various phases, the scalar field is always proportional to 
the dominant component, just as in the uncoupled model 
of ref. |(| . The three eras are clearly visible: first, the 
energy density is dominated by the radiation, with a con- 
stant contribution from the scalar field and a vanishing 
one from dark matter; then, the baryons overtake the 
radiation, and the scalar field scale accordingly; finally, 
the system falls on the final stationary accelerated at- 
tractor, where dark matter and dark energy share the 
energy density and the baryons decay away. The two 
parameters (3 and /i are uniquely fixed by the observed 
amount of fl c and by the present acceleration parameter 
(or equivalently by w e ff )■ For instance, f2 c o = 0.30 and 
W e ff = 0.33 gives fj, = 8 and P = 16, values which have 
been used in Fig. 2 and Fig. 3. With this value of /x 
we have during radiation fl^ = 6/fi 2 ~ 0.09, compatible 
with the nucleosynthesis constraints (see e.g. J|, p7|). 
To be more conservative, values of [i bigger than 11.5 
would satisfy the requirement of having < 0.045 dur- 
ing nucleosynthesis as suggested in |2j| but the situation 
would be qualitatively similar. Once (3 , [i and the present 
baryon and radiation abundances are fixed, the model is 
completely determined, and the ratio of dark matter to 
dark energy is independent of the initial conditions. 

The radiation equivalence occurs at a redshift z eq given 
by (1 + z eq ) = flbo/^l-yo — 500 for realistic values. The 
dark equivalence redshift Zdark can be found equating 
the baryon density and the dark energy density. From 
the conservation laws 



PB 



,Pc~P<j> 



a 



putting r = (3 III and approximating 
for P, ii 3> 1 ) it turns out that 



Zda 



n M (i + r) 



(14) 
(valid 



(15) 



For r ~ 2, we obtain Zdark — 5. 

The three main observations we compare our model 
to, nucleosynthesis, structure formation and present ac- 
celeration, are produced in turn during the three eras. 
The background trajectory discussed above passes the 
nucleosynthesis constraint, yields the observed accelera- 
tion and explains the cosmic coincidence. However, as 
we show next, it gives an exceedingly large integrated 
Sachs- Wolfe effect. 



III. PERTURBATIONS 

Close to the critical points bb and b c the perturbations 
in the cold dark matter component (S c ) and in the bary- 
onic one (5b) grow as a power of the scale factor, that 
is 5 C = 5b I b = a m (see ref. |Q). In this expression the 
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Table I: Critical points. 
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Table II: Properties of the critical points. 
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Figure 1: Parameter space. Each region is labelled by the 
point that is a global attractor there. Within the gray region 
the attractor is accelerated. 



Figure 2: Trend of the radiation (dashed green line), dark en- 
ergy (thick black line), dark matter (thin red line) and baryon 
(dotted blue line) density fractions, for /i = 8 and /3 = 16 
(here and in Figs. 3 and 6 the abscissa is log 10 a.) 



baryon bias b and the growth exponent m depend only on 
the parameters /_t and 0. Considering only the dominating 
growing modes, during the baryonic phase the perturba- 
tions evolve asymptotically with the law (see for instance 
ref. §) 



7711 



1 



25- 



108 



(16) 



while in the last plateau the common growth exponent is 
1 



7772 



4(/3 + /x) 



-10/9-/1 + A] 



where 



A 2 = -108 + 44/9/t + 32/3 3 /i + 25/i 2 + 1 (32/i 2 - 44) 
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Figure 3: The effective equation of state to e // = 1 + ptot/ ptot 
for the same parameters as in the previous figure. Below the 
dashed line the expansion is accelerated. 

In ref. |}4| it is also reported how restrictions on the 
baryon bias would result in further constraints on the 
model but these will not be used here since experimen- 
tal bias determinations still remain rather uncertain. The 
constraint provided by nucleosynthesis can surely be con- 
sidered on firmer grounds. 

The nucleosynthesis constraint fi > 7 (so that 0^(1 
MeV) < 0.1) together with the limitation 0.6 < fl^o < 
0.8 implies a value of (3 comprised between 9.8 and 27.3. 
For this range of values the growth exponent in the last 
era is found to lie between 7.4 and 15.3 respectively. In 
Fig. 4 numerical evolutions of 6 C and 6b using the full set 
of equations are shown in the case \x = 7 and j3 = 9.8 for a 
fluctuation of wavelength 10 Mpc hT 1 : the fast growth 
during the final stage shows up clearly. With such a 
conspicuous growth one expects a very large (late) ISW 
effect on the CMB, which in fact appears in the numeri- 
cal integrations of the model of Fig. 5, produced using a 
version of CMBFAST modified for the dark energy. The 
angular power spectrum is forced, by the normalization 
procedure, to be highly suppressed at the lowest angular 
scales respect to the observed values. This determines the 
failure of the model investigated here. Even neglecting 
the nucleosynthesis constraint, the existence of a radia- 
tion era requires fi > V6, a value that induces again an 
unacceptably large ISW effect. 

IV. CONCLUSIONS 

Our coupled dark energy model provides a cosmolog- 
ical scenario with quite unusual features. The standard 
sequence of a dark matter era followed by a cosmolog- 
ical constant era, which forces to put the present uni- 
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Figure 4: The evolution of a 10 Mpc//i perturbation for (3 = 
9.8 and /i = 7. The baryon fluctuations 5b are represented by 
the dotted (blue) line, the CDM 5 C by the continuous (red) 
line. 
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Figure 5: CMB power spectra Ce for two set of parameters: 
(3 = 4,/i = 3.5 (top curve) and (3 = 3.3,/i = 4.2 (bottom 
curve), compared with observational data from Boomerang 
|fo| and COBE The other input values for CMBFAST 

are h = 0.8, fi 6 = 0.04, fi c = 0.3, n = 1. The strong ISW effect 
at small multipoles is evident. Larger values of fi, as needed 
from the nucleosynthesis constraint, enhance the problem. 



verse on a unlikely transient, is here replaced by a bary- 
onic era and a stationary dark era in which dark energy 
and dark matter share a constant fraction of the total 
density. Contrary to almost all models published so far, 
the present universe can be seen as already on the final 
global attractor (except that, luckily, some baryons are 
still around). The two new dimensionless constants in- 
troduced in our scenario, [3 and fi, are determined by the 
present dark matter energy density and by the present 
acceleration, and can be of order unity. All cosmologies 
with an accelerated stationary global attractor reduce 
asymptotically to the model discussed in this paper. 

In this model the coincidence problem is immediately 
solved by setting (3 and fi to the same order of magnitude. 
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Regardless of the initial condition, the universe evolves 
to a stationary state with f^/f2 c =const. and of order 
unity. 

Moreover, this model explains also why the accelerated 
epoch occurs just before the present or, equivalently, why 
there are far less baryons than CDM. The reason is pro- 
vided by Eq. ([i"5|): fixing r = of order unity and Qbo 
of the order of a few per cent, we have Zdark near unity, 
regardless of the initial conditions. That is, the fact that 
we observe a relatively small quantity of baryons around 
implies that the accelerated epoch is recent. Much more 
or much less baryons would push the beginning of the 
accelerated epoch far in the future or in the past. 

The problem of the near coincidence between the ra- 
diation equivalence and the dark equivalence can be 
rephrased as why z eq and Zdark are relatively close to 



each other. The answer is that is the end of the radia- 
tion era that triggers the onset of the baryon era, which 
in turn lasts for a relatively short time because the sys- 
tem is heading toward the global attractor represented 
by the dark era. 

Despite these positive features, the model as it stands 
cannot explain our universe, since baryons and CDM 
fluctuations grow excessively during the last accelerated 
phase if it is considered that a standard nucleosynthe- 
sis has taken place earlier. This causes a CMB angular 
power spectrum suppressed at the lowest angular scales 
respect to the observational results. Thus we are forced 
to conclude that, for the universe to fall on the station- 
ary attractor, a non-linear modulation in the coupling 
(as in ref. |2^]), and/or a potential that reduces only 
asymptotically to a pure exponential, is needed. 
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